A Variant Approach to the Overlap Action 
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I describe an implementation of the overlap action, which is built from an action which is itself 
an approximate overlap action. It appears to be about a factor of 15-20 less expensive to use, than 
the usual overlap action with the Wilson fermion action as its kernel. Ingredients include a fat 
link to suppress coupling to dislocations and a free field action with a spectrum which resembles an 
overlap; much of the gain comes from the use of eigenmodes of the approximate action to begin the 
overlap calculation. As a physics example, I compute the quark condensate in finite volume in the 
quenched approximation. 



I. INTRODUCTION 
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' There has been a great deal of interest recently in lattice fermion actions which implement an exact chiral symmetry 
without doubling jl] . The one explicit realization of such an action is the overlap action of Neuberger . It obeys 
I ' the simplest version of the Ginsparg- Wilson (G-W) relation. All published studies of the overlap action in four 
^ , dimensions use the Wilson fermion action as their starting point. These actions are apparently very expensive to 

simulate. However, it is a common expectation that a better starting action would simplify the evaluation of an overlap 
action, if many fewer steps of iteration would compensate for the increased cost per step of the more complicated 
action. That expectation is realized by an action (or family of actions) I describe here. With it, some physics 
' calculations using an overlap action in quenched approximations can be done without recourse to supercomputers. 
, Ingredients in my scheme include a fat link to suppress coupling to dislocations and a free field action with a spectrum 
' which resembles an overlap; much of the gain comes from the use of eigenmodes of the approximate action to begin 
, the calculation of overlap eigenmodes. 
^ ' These techniques are obviously inspired by the fixed point action program for constructing classically perfect fermion 
I ' actions ]l3[ . The best action I have found is, however, not, as far as I know, a fixed point action of any renormalization 
group transformation. 

f-H \ Versions of these ideas have been presented many times in two dimensional models [p^ 



After setting some conventions in Sec. O, I will describe the candidate actions and their tests in Sec. 111, I will 



then present in Sec. IV a calculation of the quark condensate in finite volume with one of the new actions, along the 



lines of the calculation of Ref. f[ll| 



II. NOTATION AND CONVENTIONS 



I will call a generic lattice Dirac operator d and the overlap Dirac operator D. The eigenvalues of the simplest 
implementation of a G-W action lie on a circle of radius tq and the (massless) overlap Dirac operator is 

i?(0) = r„(l + ^) (1) 

where z = d{—ro)/ro = {d~ ro)/ro and d{m) ^ d + m is the massive Dirac operator for mass m (i.e. rp is equivalent 
to a negative mass term.) The overall multiplicative factor of tq is a useful convention; when the Dirac operator d is 
thought of as "small" and Eqn. |l] is expanded for small d, D ~ d. Apart from this overall factor of tq, my conventions 
are those of Ref. . 

The Hermitian Dirac operator for mass m is defined as h{in) — "f^d{m) and the overlap Hermitian Dirac operator 
is denoted as H{m). Specifically, 

H{Q) = ro{i^ + e{h{-ro))) (2) 
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where e{x) is the step function, e = —1 if a; < 0, e = 1 if x > 0. I will refer to the argument of the step function as 
the "kernel" of the overlap. 

It is convenient to define the massive overlap Dirac operator as 

771 

D(m) = (1 )D(0)+Tn (3) 

2ro 

and then the squared massive Hermitian Dirac operator is 

2 

H{m)^ = + (1 - ^)i/(0)2. (4) 

The zero eigenvalue eigenmodes of H{0) are chiral, with {4>\'y5\4>) = ±1, and the nonzero eigenvalue eigenmodes of 
H{0) come in pairs of equal and opposite eigenvalues. For these modes, (0a ItsI^a) = {4'>.\H{0)\(j)x} /{^ro) = A/(2ro). 

In a background gauge field carrying a topological charge Q, D{0) (and H{0)) will have Q pairs of real eigenmodes 
with eigenvalues and 2rQ. In computing propagators (for example for (^'ipip'^), it is convenient to clip out the 
eigenmode with real eigenvalue 2ro, and to define the subtracted propagator as 

Dim)-' = ^^^[D{m)-' ~ (5) 



III. THE OVERLAP WITH A BETTER KERNEL 



A. Ingredients 

A good dor h for use as a kernel in Eq. [l] or || should already "look like" D or H. This means that the eigenvalues of 
its eigenmodes should lie (approximately) on a circle and its low lying eigenmodes should be approximately chiral-in 
a spectroscopy calculation, the additive renormalization of the bare quark mass, as measured, for example, through 
the variation of the pion mass with bare quark mass, should be small. Most implementations of the overlap action in 
the literature use the Wilson action, which does not satisfy either of these criteria: the eigenvalues of the free Wilson 
action sit like beads on a string around a set of four circular arcs (with real parts of their eigenvalues ranging from 
to 8), and the additive mass renormalization of the interacting theory is Aamq ~ 1 for simulations at lattice spacings 
near 0.15 fm. Use of the clover action instead of the Wilson action improves the chiral properties but does nothing 
for its (free field) eigenmode spectrum. 

Improving the action involves two ingredients. 

First, to improve the chiral properties of the action, thin links are replaced by fat links. Actions with fat links are 
already quite chiral as shown by their small mass renormalization and Za — 1 ||l5| , |l6|] . In this work I have studied fat 
link actions with two blockings: The first uses APE- blocking | |T7[ |: 

+ v!,''-'\x - v)^vl.''-'\x - J^)v!,''-^\x - V + A)), (6) 

with vIP^ [x) projected back onto 5[/(3) after each step, and {n) — U^{n) the original link variable. I have mostly 
studied a large amount of fattening, a = 0.45, N = 10, but have looked at the smaller value a = 0.25, N = 7. The 
second blocking is a "hypercubic" blocking devised by A. Hasenfratz jl8|. The fat links of this action are confined to 
a hypercube, so this action is more local than the APE-blocked actions. The mean plaquette at /3 = 5.9 for (0.45,10) 
APE blocking is about 2.98, for (0.25,7) APE blocking it is 2.88, and for the hypercubic blocking, 2.84. 

I believe that what is important here is that the fat gauge links decouple the fermions from short distance structure 
in the gauge field, not that the links are fattened in a particular way. 

The eigenvalues of a GW action lie on a circle. I determine the best d by taking a free field test action and 
varying its parameterization to optimize its eigenvalue spectrum (in the least-squares sense) for circularity, for some 
rp. It was convenient to let the parameter vq also be a free parameter. The simplest extension of a nearest neighbor 
action which can have its eigenvalues lying approximately on a circle is "planar:" it has scalar and vector couplings 
S = J2x r '4'ix){^{r)+i"f^Pf^{r))'il!{x + r) for r connecting nearest neighbors (r — ±fl; A = Ai = —0.170, p'j}^ — —0.177) 



2 



— * (2.) (2') 

and diagonal neighbors (r = ±/i ± i>, 7^ /x; A = A2 = —0.061, —pi = —0.0538. The constraint A(r = 0) = 

— 8A1 — 24A2 enforces masslessness on the spectrum, and — 1 = 2p^i}^ + \2p^^^ normahzes the action to —Tpi^p^d^ip in 
the naive continuum hmit. The corresponding value of tq is 1.6. A plot of the free field eigenvalues on a finite lattice 
is shown in Fig. |^. Presumably the places where its eigenvalues are purely real (but not zero) correspond in the 
continuum limit to theories with various numbers of massless free fermions, exactly as for the Wilson action. 




Re A 

FIG. 1. Eigenmode spectrum for the free field planar action on a finite lattice. Only the positive imaginary part modes are 
shown. 

The action also includes a clover term, with its coefficient set to the tree-level value appropriate to this action of 
Csw = 1.029941. Thus d has no 0{a) artifacts. All the links, including the ones in the clover term, are replaced by 
fat links. The gauge connections to the diagonal neighbors are the average of the two shortest path connections. The 
cost of the action is about 6.5 times that of the usual clover action. 

While it is probably not germane to the present discussion, the dispersion relation of this action is improved 
compared to that of clover or Wilson actions. 

No claims are made for uniqueness-or even optimality. The optimizations do not have to be done with enormous 
precision since the overlap formula itself deforms almost any kernel action into a chiral action. The amount of 
fattening is also a free parameter, and there will almost certainly be tradeoffs between ease of implementing the 
overlap and the desire that the action show good scaling behavior. At large levels of fattening, there appears to be 
little renormalization of the parameters of a free action, and so little tuning is required to produce actions which 
behave well in simulations. Thus, the overlap parameter ro will be kept at its free field value of 1.6. 
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B. Testing actions 



In what follows I will focus on comparisons of the overlap action with the Wilson action as a kernel, the "Wilson 
overlap" and with the fat link planar action as a kernel, or "planar overlap." I have looked at six actions in all. Most 
tests involve the Wilson overlap or the (0.45,10) APE-blocked planar overlap. I have briefly investigated an overlap 
action with the (0.45,10) APE-blocked fat link clover action as a kernel. I have also investigated planar overlaps with 
(0.45,10) APE-blocked links and with hypercube-blocked links. Finally, I tested the "Gaussian" action, which has 
been proposed as a candidate approximate FP action [ p^ . 

The step function is approximated by either the polar formula introduced by Neuberger 

1 ^ 1 

eiz)^e^iz) = z-Y.^;^ (7) 

{cj = cos'^(7r(j -I- l/2)/{2N)), sj = sin^(7r(?' + l/2)/(2A^))), or by a fourteenth order Remes polynomial, following the 
work of Edwards, Heller and Narayanan [fj. In practice, the polar formula works very robustly for the planar overlap, 
when I rescale h{—ro) to h{~ro)/ro in Eq. ^ and take = 6 to 10. The Wilson overlap is much more delicate and I 
have used the Remes approximation exclusively for it. I rescaled the operator h{—ro) by a factor of 1/2.5 to map its 
eigenvalues into a range where the Remes algorithm has small errors. The inverses of the terms in the sum of Eq. 
are found using a multi-mass conjugate gradient routine p9| . 

Let's begin by looking at the actions. The range of the action is computed using a variation of the calculation of 
Hernandez, Jansen, and Liischer [0 : Shown in Figs. || and || is a comparison of Wilson and planar overlap •\/|Z3(0)xp 

for a delta- function source x at the origin, as a function of distance r — ^J^, for the free field action as well as on a set 

oi P = 5.9 8^ configurations. The two actions show similar exponential falloff with distance. Of course, these pictures 
do not show anything about the locality of the gauge connections in the actions, just the fermions. In perturbation 
theory |Q the fat link action has a cloud of glue of size i/aA/S convoluted over every fermion offset. 
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FIG. 3. Range of the (0.45,10) APE-blocked planar overlap action, in free field theory and on a set of /3 = 5.9 8'* lattices. 

It might be relevant that the plaquette for the fat hnk action has a value Tr(l — Up)/3) ~ 1 — 2.97/3 = 0.01, which 
should be sufficient for the argument of Hernandez, Jansen, and Liischer to insure the locality of the overlap action. 

All my tests of the overlap begin with finding the eigenmodes of H with the smallest eigenmodes (actually the 
eigenmodes of H diagonalized in a basis which is composed of the smallest eigenmodes of H^). Low lying eigenmodes 
of h{m) and H{0) are found using an adaption of a conjugate gradient algorithm of Bunk et. al. and Kalkreuter and 
Simma pT[ |; I modified a code originally written for staggered fermions by K. Orginos. 

I have implemented most of the standard tricks for efficient evaluation of eigenmodes oi H . I compute a set of 
the A^o low lying eigenmodes of h{—rQ) (typically A^o = 10 — 20), and project them out during the calculation of the 
operator e{h{—ro)))x- 

A first sign that the fat link will make the overlap better behaved comes from looking at the eigenmodes of the 
kernel function /i(— tq). Fig. || shows a histogram of the ten smallest eigenmodes of ft.(— tq) reconstructed from the ten 
smallest eigenmodes of /i(— rg)^ on a set of ten S** /3 = 5.9 configurations. The fat link actions have many fewer small 
eigenmodes. Eigenmodes near zero are associated with dislocations, gauge configurations on which topological objects 
are about to disappear. It is difficult for approximations to the step function to process these modes. The absence 
of these modes in the spectrum of /i(— tq) is why the polar form of e^r can perform well, even for small N. There is 
a lot of discussion in the literature (see for example Ref. [^^ ) about improving overlap or domain wall fermions by 
using gauge actions which have fewer dislocations than the Wilson gauge action, but a fermion action which could 
not see dislocation would work just as well. A fat link action does just that. While the fat link clover action is as 
much as an improvement as the planar action, it has a much higher conditioning number than the planar action. This 
is reflected in the amount of work needed to find the eigenvalues (and will also affect the number of inner conjugate 
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gradient steps needed to evaluate the step function): the average number of Rayleigh iterations needed to collect the 
ten lowest eigenmodes is 3600 for the Wilson action, 5100 for the fat link clover action, and 1900 for the (fat link) 
planar action. 
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for the ten smallest eigenvalues of h{—ro) for (a) the Wilson action (ro = 1.65), (b) the fat 



and the fat link clover action (ro — 1.0), on 10 8 (3 = 5.9 configurations. 



Because [i?(0)^,75] = 0, eigenmodes of H{0)^ can be simultaneously eigenmodes of 75. All eigenmodes are con- 
structed in a basis of chiral eigenfunctions. In that basis, eigenmodes of iJ(0)^ are expected to be doubly degenerate. 
This degeneracy can be used to test algorithms. A low order polar approximation to the step function does a poor 
job of producing degenerate pairs of eigenmodes of H{0)^ for the Wilson overlap and at this level of test it must be 
discarded. 

The calculation of the condensate in Sec. IV requires the quark propagator D(jn)^^, with the contribution from 
chiral modes removed. 1 do this by taking the set of eigenmodes of h and finding the two lowest eigenmodes 77(0)^ 
in a basis containing one state of each chirality. The chirality of the highest state identifies the chiral sector which 
has no zero modes. I then perform the inversion of = D'' D in this sector. This calculation is also accelerated by 
finding and projecting out some large number (0(10)) of eigenmodes of H{0)^ during this inversion. 

A trick which I have not noticed to have been emphasized in the literature is actually the source of much of the gain 
in efficiency of the planar overlap. This is to begin the calculation of eigenmodes of H by first finding a set of low lying 
eigenmodes of h{m) for some potentially useful m, and beginning the calculation for H using the eigenmodes of h 
(rather than, say, beginning with a set of random vectors). The combination of the small number of small eigenmodes 
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of h{—ro) for the fat link action and the use of good trial functions makes the planar overlap quite efficient compared 
to the Wilson overlap. 

Most of my tests are on a set oi8* (3 = 5.9 lattices. Let us consider a set of examples which illustrate the differences. 
In the planar overlap I take the A'^ = 10 polar approximation to the step function, scale the argument by — 1.6, 
and project out A^'o = 10 eigenmodes of h{—ro) in the step function. The Wilson overlap uses the fourteenth-order 
Remes approximation, scales its argument by 1/2.5, and projects 20 eigenvalues of h{~ro). In both cases I wish to 
find the lowest two eigenmodes of H{0)^, regardless of chirality, and I begin by finding and utilizing the four smallest 
eigenmodes of h{0). I investigate a configuration which happens to have topological charge Q = 1. The planar overlap 
calculation needs 109 Rayleigh iterations and 2529 inner conjugate gradient steps to find the two lowest modes (at a 
cost of 6.5 times the equivalent Wilson action step). The Wilson overlap needs 1220 Rayleigh iterations, and 86,723 
inner conjugate gradient steps, to reach the same level of accuracy (about 10~^ for the eigenvalues). Had we wanted 
the ten smallest eigenmodes of H{0)^, using 20 trial eigenmodes of h{0)^, the cost would be 306 Rayleigh iterations 
and 8189 inner conjugate gradient steps for the planar overlap, or about 26 inner steps per Rayleigh iteration, while 
the Wilson overlap uses 5820 Rayleigh iterations and a million inner conjugate gradient steps (about 170 inner steps 
per iteration) . This is a savings in computer time of about a factor of eighteen for the planar overlap compared to the 
Wilson overlap. On other configurations, the number of inner CG's per Rayleigh step for the planar overlap ranges 
between 20 and 30, and the number of Rayleigh iterations ranges from 300 to 500. 

The difference in the number of inner conjugate gradient steps is due to the fact that the Wilson h{~ra) has many 
more small eigenvalues than the fat link planar action. The difference in the number of Rayleigh iterations arises 
because the initial trial vectors, eigenmodes of h{0), are close to being eigenmodes of -ff(O). I illustrate this with a 
scatter plot of the change in an eigenmode, plotted as a function of the value of the overlap eigenmode, in Fig. ||. 
A graph of (75) vs. eigenmode for this action, at various stages of the overlap calculation, is shown in Fig. ^. It is 
also useful to look at the scatter of chirality vs. eigenvalue for several trial actions which might be candidates for an 
overlap kernel. This is shown in Fig. ^. I show (75) vs. eigenmode for the Wilson action, non-perturbative thin link 
clover action (with Csw — 1-85), fat link planar action, and fat link Gaussian action. From this picture, one would 
suspect that the thin link actions would not provide good trial wave functions for the overlap. 
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produced by diagonalizing H{0) in the lowest 20 eigenvalue eigenmodes of h{0)^, as a function of overlap eigenvalue. 
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FIG. 6. {75) vs. eigenvalue for the planar overlap: (a) 20 smallest eigenvalues of /i(0), (b) 20 smallest eigenvalues of -ff(O) 
from a diagonalization of the basis of (a), (c) 20 smallest eigenvalues of H{0) from a full calculation, all on a set of ten /3 = 5.9 
8'* lattices. 
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FIG. 7. (75) vs. eigenvalue for several h{0)'s. (a) thin link Wilson action, (b) non-perturbative thin link clover action, (c) 
Gaussian action, (d) planar action. 

I have tested another fat link action, a planar action with 7 APE-blocking steps and a = 0.25. I kept rp = 1.6 since 
tests of the fat link clover action with this level of fattening at /3 = 5.9 also showed little additive mass renormalization. 
The number of small eigenvalues of h{—ro) increases slightly compared to the (0.45,10) case, and the number of inner 
CG steps grows, correspondingly, to about about 30. But the eigenmodes of h{0) do not seem to be as good a 
trial basis as they were for the fatter link action, and the number of Rayleigh iterations grows from about 10,000 to 
12,000-15,000 for the same calculation as was done above. A set of pictures showing relevant information for this 
action is shown in Fig. |[ This action does not see some of the instantons that the (0.45,10) fattening saw. That is 
reflected in the scatter of the eigenvalues. 
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FIG. 8. (75) vs. eigenvalue for the planar action and its overlap with (7,0.25) fattening: (a) 20 smallest eigenvalues of 
h{0), (b) 20 smallest eigenvalues of H{0) from a diagonalization of the basis of (a), (c) histogram of ten smallest eigenmodes of 
h{—ro), on /3 = 5.9 8^ configurations. 

One might hope that a better trial function would help. A simple way of changing the trial wave function is to vary 
the mass m in h{m). Trials at m = —0.05 and —0.1 did not produce any dramatic changes. This is an obvious place 
for future work. 

The planar action with hypercubic fat links shows nearly identical behavior. Both these actions need about fifty 
per cent more inner CG steps than the (0.45,10) planar action but are still a large improvement over the Wilson 
overlap-and one could still run them on work stations. They might be useful in contexts where a very fat link might 
have undesirable properties. 

Finally, the Gaussian action with (0.45,10) APE blocking requires about half as many inner CG calls as the (0.45,10) 
planar action, typically 3000-5000. However, this action requires about 2.5 times as much storage and 2.5 times as 
much CPU time per inner CG step as the planar overlap action, since its couplings span a hypercube (80 neighbors 
rather than 24) . Thus it does not seem to be competitive with the planar action without clever coding |23| . 

The (0.45,10) fat clover action, with Csw = Ij was also briefly investigated as a kernel for the overlap. It seemed 
to need 100,000-300,000 inner CG steps in the fiducial calculation, in 2500 to 8000 Rayleigh iterations, beginning 
from m — trial wave functions. (There is an extreme variation in these numbers from lattice to lattice tested.) 
Eigenmodes of the fat link clover action are apparently quite different from those of the corresponding overlap action. 
The number of inner CG steps reflects the greater conditioning number of /i(— tq). 

These tests are certainly incomplete. It may be that I have simply written a very inefficient Wilson overlap, though 
most of the code for the two methods is common. My Wilson overlap is so expensive that it cannot be tuned, whereas 
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it is easy to test variations of the planar overlap. The "design philosophy" of beginning with an action which is 
close to an overlap action certainly produces actions which are inexpensive enough to run on small computers (work 
stations in my case). 

IV. THE QUARK CONDENSATE IN FINITE VOLUME 

In order to do a little physics in what is otherwise a paper about technique, I present a calculation of the chiral 
condensate in the quenched approximation with the (0.45,10) fattened planar overlap action. The method is (almost) 
exactly that of the pioneering calculation of Hernandez and Jansen and Lellouch |ll| : one computes the condensate 
in background gauge field configurations of fixed topology labeled by winding number v . The condensate is 

= raY?V{I^{mY.V)K^{mLV) + /^+i(mI]T/)if^_i(mST/)) + (8) 

where m is the quark mass, V is the volume, and E is the infinite volume condensate. lv{_z) and K^iz) are modified 
Bessel functions. In practice, one computes E "without topology," by working in the chiral sector which has no zero 
eigenmodes (and doubling the result). Then one measures E^ — vjimV^. In practice, the lattice number needs an 
additive renormalization: in zero mass, this is removed by the replacement of by D^^ as in Eq. I There can 
also be contribution which vary with the quark mass, so the lattice number will need to be fit to 

- E ^--'^ = - ^lij^y)) + Cm + ... (9) 

X 

One small difference between my calculation and that of Ref. (ll| is in the normalization of the propagator: they do 
not have the prefactor 1/(1 — ^) of Eq. ^. It is hard to believe that the condensate should not be an odd function of 
the symmetry breaking term (the quark mass), and the data does not show anything but a simple linear dependence 
on the mass. 

Finally, E is scheme-dependent, and an overall lattice-to-continuum regulator needs to be computed. 

Warned by the results of Ref. ||ll[| , I restricted my calculation to the v = ±1 sector. I generated a set of 40 lattices 
at each of three volumes S**, 10^, and 12'* at /? = 5.9. I filtered them to find candidate v = ±1 configurations using a 
pure gauge measurement of topological charge which I had previously calibrated against a set of 10 8"* lattices. 
I checked that these lattices in fact had v — ±1 during the calculation of the condensate; only one lattice failed this 
test. This obviously leaves out configurations which the fermion observable identifies as carrying topological charge, 
but the gauge observable does not, but the alternative is to process every configuration through the overlap program. 
I was left with a set of 10 S**, 13 10'*, and 9 12^^ lattices on which I computed I computed propagators at five 

quark masses, am = 0.001, 0.002, 0.005, 0.01, and 0.02 (again using a multi-mass Conjugate Gradient solver). I used 
twelve random sources per lattice (though I blocked data together before averaging). This needed about 40, 70, and 
125 steps at the three volumes, at about 20-30 inner CG's per outer step, to reach a fractional squared residue of 
10~*^. As a check, I also calculated E from the GMOR relation — ^/"tn; as expected in the overlap it 

reproduced the direct calculation of E/m to within 10^^ (see Ref. 0). 

My lattice results are shown in Fig. |[ 



12 



0.55 



0.50 



0.45 



0.40 



0.35 



0.30 



X 



J I I L 












i 


> i 




¥ — , 




= E 


3 [ 


] [ 


3 E 














3 E 


J- 



J I I I I I L 



1x10"^ 2x10"^ 



5x10"^ 1x10"^ 2x10"^ 



m 

q 

FIG. 9. The quark condensate in finite volumes (divided by the quark mass). From the bottom, the volumes are 8*, 10*, 
and 12*, The curves are fits to Eq. |l^. 



The data at each volume are quite correlated. A single-elimination jackknife fit to 

(-0^;)^^^ = mi:'^V{Ii{mY.V)Ki{mJ:V) + l2imT.V)KoimT.V)) + Cm (10) 

gives (re-inserting the lattice spacing) Sa"^ — 0.00394(16], C — 0.304(8). The lattice number is about two standard 
deviations higher than the number reported by Ref. |1^ at /3 = 5.85 of Sa'^ = .0032(4). We can sharpen this 
disagreement by trading the lattice spacing a for the Sommer radius tq, using the interpolating formula of Ref. [ p6| : 
Erg = 0.215(27) for Ref. |l|], Erg = 0.354(14) here. 

A real comparison requires computing the Z factor. I have not done that yet, but I can make a heuristic attempt 
at the calculation, by exploiting the fact that perturbation theory for fat link actions becomes sim ple in the limit of 
large fattening. This argument is implicit in the discussion of fat link perturbation theory in Ref. pG|] . Since it falls 
outside the main thrust of the paper, I relegate the discussion to an Appendix. 

Taking rg to be 0.5 fm, ^^j^ifJ. = 2 GeV) = 0.0244(10X37) GeV^. The two errors are from the lattice fit and from 
an assumed lattice spacing uncertainty of five per cent [E^ . 

This calculation produces a number which disagrees badly with a calculation of the condensate using clover fermions 
and the GMOR relation, ^jjgifJ- = 2 GeV) = 0.0147(8) (16) (12) GeV^ It is done at smaller lattice spacing and has 
completely different systematic uncertainties. Oddly enough, however, my result agrees well with an earlier analysis 
of data by Gupta and Bhattacharya: Ep-g(/i = 2 GeV) = 0.024(2) (2) GeV^ (statistical and lattice spacing 
uncertainties). 
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Of course, all the potential weaknesses of this calculation, reported by Ref. |ll|| apply here, too: the lattices are 
small, and the lattice spacing is large. And it is a quenched calculation. 

All of the computations of S were done over about a month of running on a set of four or five 450 Mhz Pentium II 
work stations. 



V. CONCLUSIONS 



None of the results shown here are particularly surprising, but that does not mean that they might be totally devoid 
of interest. A good approximate overlap action is easier to convert into an exact overlap, than any random action 
like the Wilson action. Crucial ingredients seem to be some kind of gauge connection which suppresses dislocations, 
a free action which "resembles" a free field overlap action, and the use of as much information from the trial h{m) as 
possible to begin the calculation of the overlap action. One should be able to do better. 
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APPENDIX A: Z-FACTORS FOR VERY FAT ACTIONS 

The evaluation of the Z-factor is straightforward. I will use many of the techniques of Ref. |2^, but collect them here 
for completeness. We expect Yj{fj,)j^ = Z{aiJL)Yi(a)iatt with Zs = 1/2'™ and is the quark mass renormalization 
factor. In perturbation theory, Z„i = \ + bz, and 

z = A;ott - Acont = 6 ln(l/ (yua)) + x, (Al) 

where b = as(g*)/(37r). I choose to use as = a-jjg{q*). I will use the plaquette to define a coupling Q!y(3.41/a) 
(=0.16054 at /3 = 5.9), which will then be converted using two-loop perturbation theory to ajjg{q*), to give = 
\/a)-j-fg. I will then run this result to fi = 2 GeV, using the (inverse of the) two-loop running formula for the MS 
quark mass. 

z is the difference between a lattice integral 



k 



^latt = I hattFN] (A2) 

and a continuum integral 



Acont — / Icont (-^3) 
Jk,c 

I have factored the lattice integrand into a thin link piece Iiatt and the form factor of the fat quark-gluon vertices F^] for 
APE-smeared links Fn = {l-af/Qf^ where f = I],,(4/«^) sin2(fc^a/2). I am using the notation = / d'*fc/(27r)** 
over the hypercube, — vr/a < < tt/o; /j. ^ = / (i^fc/(27r)^. Introducing a giuon mass A to regularize its IR divergence, 

^cont = 5/2 — 31n(A^/i^) (in MS). Adding and subtracting continuum-like terms to isolate the divergence in the 
lattice integral, we write 

Aiatt - J^^illatt - 3 p )I<N --iJk p (1 - ^n) + 3 Jj^ p 

J -3D + 3ln^. (A4) 

The first two terms (J and D) are IR and UV finite. Now the point is that is a continuum-like integral, and all 
the dependence of Aiatt on the lattice action is in the first term. But if the fattening is large, the form factor F^, 
is nonzero only at tiny k, where it is unity. At that point, for all practical purposes, Iiatt — 3/k^, and so we expect 
that J will be small. All the dependence of Z„j on the lattice action, other than its form factor, vanishes when we set 
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J = in Eq. A4. q* (defined according to the prescription of Lepage and Mackenzie |29|) is calculated by a similar 
procedure. 

While I can't check this approximation for the overlap action, at (0.45,10) fattening, a complete calculation for 
clover fermions gives x ~ —5.58 and q* ~ 1.07, while the J = result is a; = —5.44 and q* — 1.10. A similar argument 
would predict that finite rcnormalization factors, such as the vector and axial current renormalizations, are unity, 
which is a good approximation to what is seen at large fattening. Of course, for small fattening, this approximation 
fails badly and a better (non-perturbative?) calculation is necessary. (For example, the (0.25,7) APE-blocked clover 
action has x = —0.975, q* = 2.88; the approximate result is a: = —2.49, q* = 1.59.) 
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